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Abstract 

The asserted dominant role of the kidneys in the chronic regulation of blood 
pressure and in the etiology of hypertension has been debated since the 
1970s. At the center of the theory is the observation that the acute relationships 
between arterial pressure and urine production — the acute pressure-diuresis 
and pressure-natriuresis curves — physiologically adapt to perturbations in 
pressure and/or changes in the rate of salt and volume intake. These 
adaptations, modulated by various interacting neurohumoral mechanisms, 
result in chronic relationships between water and salt excretion and pressure 
that are much steeper than the acute relationships. While the view that renal 
function is the dominant controller of arterial pressure has been supported by 
computer models of the cardiovascular system known as the "Guyton-Coleman 
model", no unambiguous description of a computer model capturing chronic 
adaptation of acute renal function in blood pressure control has been 
presented. Here, such a model is developed with the goals of: 1 . representing 
the relevant mechanisms in an identifiable mathematical model; 2. identifying 
model parameters using appropriate data; 3. validating model predictions in 
comparison to data; and 4. probing hypotheses regarding the long-term control 
of arterial pressure and the etiology of primary hypertension. The developed 
model reveals: long-term control of arterial blood pressure is primarily through 
the baroreflex arc and the renin-angiotensin system; and arterial stiffening 
provides a sufficient explanation for the etiology of primary hypertension 
associated with ageing. Furthermore, the model provides the first consistent 
explanation of the physiological response to chronic stimulation of the 
baroreflex. 
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i:/47frf=/*i Changes from Version 1 

We are grateful to the reviewers for their insightful commentaries 
and constructive criticism. 

An overriding issue in the reviews is the "phenomenological" 
nature of the model equations used in this study. Although the 
phenomenological nature of the model was not judged to be 
a shortcoming, the issue was raised repeatedly. Yet all models 
are "phenomenological" at some level. Indeed, even the models 
pointed out by the reviewers incorporate exactly the same sort of 
phenomenological representations that we have invoked. Thus, 
rather than the "phenomenological" nature of our model, the more 
specific issue may be our model's simplicity. Indeed, certain 
aspects of our model are simple compared to certain aspects 
of certain other models. That simplicity could be judged to be a 
strength or a weakness, depending on the context or application. 
Thus, it is our opinion the important questions are: Does the 
model capture the phenomena that we set out to capture? Is the 
model validated against experimental data that were not used 
to identify the model? Does the model facilitate any insights into 
physiological phenomena? In other words, we invoke the often 
repeated words of George Box, "All models are wrong, but some 
are useful". 

So, how is our model useful? The reviews have helped us a great 
deal to spell out clear answers to that question. In the revised 
text we have clarified the goals and stated the conclusions so that 
they are clearly justified by the analysis presented in the paper. In 
addition, a basic parameter sensitivity has been performed with 
results presented in Figure 1. This analysis revealed a slightly 
improved parameterization of the model, reflected in the updated 
Table 1 and Figure 6-Figure 10. Finally, we have added a panel to 
Figure 4 illustrating the overall organization of the model. 

See referee reports 



Introduction 

Theoretical analysis and observations of the control of blood vol- 
ume and salt content by the kidneys has led to the hypothesis that 
arterial pressure is determined in the long-term (over time scales 
of days or more) by the balance between the level of salt intake 
and the acute relationship between pressure and salt excretion by 
the kidneys. In fact, it is stated that the renal pressure/volume con- 
trol system adjusts arterial pressure with "infinite gain" and that the 
renal function curve and rate of salt and water intake are the "two 
sole determinants of the long" 1 . Yet while it is certain that at any 
steady level of arterial pressure the rates of salt and water intake and 
excretion are balanced, there is considerable debate over how this 
balance is achieved, and thus what are the long-term determinants 
of arterial pressure 2-5 . 

One viewpoint is that the acute relationship between pressure and salt 
excretion — the acute pressure-natriuresis mechanism — represents a 
physiological input-output relationship and that alterations to this 
mechanism underlie most (if not all) chronic changes in pressure 1,6 . 
While competing ideas identify the nervous system and its influ- 
ence on the heart and vasculature as the primary long-term con- 
trollers of arterial pressure 3 4 , the hypothesis that long-term control 
of blood pressure is achieved through the renal pressure-diuresis 
represents the dominant thinking in the field 7 . 

The most widely recognized model of long-term blood pressure 
control is the model developed over several decades by Guyton 



and colleagues. The 1972 realization of the model 6 invokes 
approximately 160 variables, and several hundred adjustable param- 
eters 8 . More recent versions of the model involve many thousands 
of variables and tens of thousands of parameters 910 . While several 
realizations of the Guyton model are or have been disseminated 
as computer programs, human readable expositions of the various 
versions of the model are lacking. Moreover, there is no published 
report of any of these models in which the governing equations 
are defined, the parameter values are provided, the data that were 
used for model identification are reported, and the model identifi- 
cation process and results are reported. Indeed, it is unlikely if not 
impossible that even the relatively simple 1972 version of the model 
was formally identified. These facts unfortunately make this seminal 
work largely impenetrable and possibly irreproducible. Most criti- 
cally, the Guyton-Coleman model is based on the assertion that arte- 
rial blood pressure is, under all circumstances, controlled primarily 
by the kidney 11 . Therefore the Guyton-Coleman model cannot be 
used to explore alternative hypotheses. 

In the Systems Approach for PHysiological Integration of Renal, 
cardiac and respiratory functions (SAPHIR) project, Thomas 
et al} aimed to develop and disseminate a core "model of human 
physiology targeting the short- and long-term regulation of blood 
pressure, body fluids and homeostasis of the major solutes" includ- 
ing "the main regulatory sensors (baro- and chemoreceptors) and 
nervous (autonomic control) and hormonal regulators (antidiuretic 
hormone, aldosterone and angiotensin)". The project goal is to, for 
the first time, provide an open- source transparent model integrating 
these systems. In 2008, the SAPHIR group published a core model 
of circulation and volume exchange across fluid compartments that 
is derived from the Guyton models, consisting of 20 equations and 
perhaps less than 100 parameters 8 . This model does not account 
for adaptation of the acute pressure-diuresis/natriuresis function to 
chronic changes in arterial pressure or volume/salt loading. Thus, 
it does not capture the concept of the chronic renal function curve, 
which is thought to be of fundamental importance to the long-term 
regulation of arterial pressure. Recently Averina et al. 4 reported a 
mathematical model of long-term control of arterial pressure that 
captures the concept of chronic adaptation of pressure-natriuresis 
and, to our knowledge, represents the only published model incor- 
porating the concept of adaptation of the acute pressure-natriuresis 
relationship in response to changes in salt/volume loading. Thus, the 
model of Averina et al. represents the model of record of regulation of 
blood pressure accounting for the chronic renal function curve. Yet, 
this model was developed to illustrate that the chronic renal function 
curve need not represent a physiological input-output mechanism, 
and to introduce an alternative description of the long-term regula- 
tion of arterial pressure and the etiology of salt- and angiotensin- 
dependent hypertension. 

Similar to the goals of the SAPHIR project, we developed a math- 
ematical model of the Guyton concept of long-term control of arte- 
rial pressure, analysed its behavior and compared its predictions 
to experimental observations. Similarly, an open- source transparent 
model integrating these systems is provided. The approach pursued 
here differs from that of Thomas and colleagues in that while the 
SAPHIR project explicitly centers on translating the components of 
the Guyton model into an explicitly defined core model, our goal is 
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to develop a model from scratch in which all model components and 
associated data are formally identified based on experimental data. 
Furthermore, given the demands of model identification and vali- 
dation, our goals are necessarily more narrowly constrained than 
those of the SAPHIR project and of the original Guyton-Coleman 
models. In our study, a model is constructed in the spirit of the 
Guyton-Coleman models, adopting a practical approach to data- 
driven phenomenological representations of physiological systems, 
rather than aiming for physical and mechanistic realism. 

Specifically, we have developed a model of the long-term control 
of arterial pressure that captures the Guyton concept of pressure- 
diuresis/natriuresis as physiological input-output relationships. 
The model is designed to meet the following criteria: 1 . The model 
accounts for the effects of the baroreflex on the circulation, the heart, 
and the kidney, and the effects of the renin-angiotensin system on 
the circulation, the heart, and the kidney; 2. The model is explicitly 
documented for definition and reproducibility: all model equations 
are reported herein and justified based on experimental data; 3. 
Parameter values are reported and justified: All parameter values are 
estimated by comparing model simulations to measured data. The 
identified model is used to probe physiological mechanisms under- 
lying the chronic renal function curve, to provide insight into how 
chronic stimulation of the baroreflex leads to chronic reductions in 
arterial pressure, and to generate hypotheses regarding the etiology 
of primary hypertension. 

Methods 

Model Components 1 : Aorta/large-artery mechanics 

Based on the simple approximation of a thin- walled cylinder, the 
strain e in the aorta is computed as a function of volume V Ao 



do 



(i) 



where V Q is a parameter representing the unstressed volume, and 
d A /d Q is the ratio of diameter to unstressed diameter. (Here e is 
defined to be equal to 1 in the unstressed state when d A = d Q ). The 
aortic pressure-volume relationship is assumed to be governed by 



-V )/c 



(2) 



where C Ao represents an acute compliance and V sAo (t) accounts for 
creep mechanics of aortic wall, simulated according to 



dV 



sA^-y™ _y 

sAo sAo 



V =y V 

v Ao / Ao y Ao> 



(3) 



where t , is the time constant of stress relaxation and C, /C is 

cAo Ao °° 

the ratio of acute to effective chronic compliance of the vessel. 
(Equation (2) and Equation (3) represent an alternative equivalent 
formulation of the standard linear model of vessel mechanics). 



The parameters of the large-artery mechanics model component are 
identified based on measurement of aortic diameter and pressure 
in dogs reported by Coleridge et al. 12 . Figure 1 plots the measured 
aortic pressure wave in the upper panel and vessel diameter as a 
function of pressure in the lower. The experiment of Coleridge et al. 
is simulated using the governing equations 



dV. A 



dt 



--V°° -V 

sAo sAo 



dV ^-C A dP *° i dV « 
dt dt dt 

d£ = 1 dV An 
dt 2{V Ao V 0 ) h 



(4) 



where the measured aortic pressure waveform (Figure 1A) is used 
to numerically approximate dP A0 (t)/dt in integrating Equation (4). 
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Figure 1. Simulated aortic mechanics. A. The aortic pressure 
time course obtained from Coleridge et al." is used as the input to 
the aortic mechanics model module, Equation (4). B. The model- 
predicted relationship between aortic pressure and diameter 
is compared to the data reported by Coleridge et al.". Model 
simulations are plotted as a solid black cure; data are plotted as 
shaded circles. 
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The parameters y Ao , C Ao , and t cA , were adjusted to match the data in 
Figure IB; the values of V^and d Q were set arbitrarily by assuming 
a cylindrical vessel of length 30 mm. All parameter values are listed 
in Table 1 . 

Model Components 2: Kinetics of baroreflex afferent firing 

The baroreceptor afferent firing rate is assumed to be governed by 
the rate of change of strain in the vessel wall. The model invokes a 
moving average strain value i{t), which is computed 



de 



-£-£ . 



(5) 



The time constant t is an adjustable parameter. The baroreceptor 
firing rate is assumed proportional to S £ = max(£-£,0) via the satu- 
rable relationship 



f BR (t)=f 0 s(t) 



S e (t) 
8 B (t)+8 0 ' 



(6) 



where s(t) represents the fraction of baroreceptor afferents that are 
in an active/permissible state, and 8 Q and/ 0 are additional adjust- 
able parameters. Equation (6) is a static nonlinearity 13 that enforces 
a saturating response. It is assumed that the baroreceptors within 
the population transition from an active to inactive state at a rate 
proportional to the firing rate, and transition to the active state at a 
constant rate: 



—=a(l- s) - 



bs- 



8 c +5 n 



(7) 



The adjustable parameters in the baroreflex afferent model (t, 
5 0 , f Q , a, and b) are identified based on measurements following 
step changes in non-pulsatile carotid sinuspressure 14 and ramps in 
in vivo pulsatile aortic pressure 12 . 

The data in Figure 2A, obtained from Chapleau et al. 14 , correspond 
to the multi-fiber response of the carotid sinus baroreceptors fol- 
lowing a step change in pressure from 40 mmHg to 93 mmHg. To 
simulate this experiment, the pressure is assumed to follow a time 
course described by 



40,f<-l 
P Ao (t) = 40 + 53- (; + l),;>-l, 
93, t>0 



(8) 



where pressure is expressed in units of mmHg and time in seconds. 
Thus it is assumed that pressure is increased by +53 mmHg over an 
interval from -1 to 0 seconds. Since the overall model does not dif- 
ferentiate between aortic and carotid baroreflex signals, data from 
the carotid pressure step are matched to the model of Equation (4)- 
Equation (7) in Figure 2A. Following the step pressure increase the 
afferent firing rate f BR rapidly increases to a maximum value before 
decaying over a time scale of several hundred seconds. This decay is 
determined in the model by the combined action of mechanical relaxa- 
tion governed by Equation (5) and inactivation governed by Equation (7). 
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Figure 2. Baroreflex model. A. The response of the baroreflex model 
to a step increase in pressure is compared to data from Chapleau 
ef a/. 13 . The model simulations are based on Equation (4)-Equation (7), 
as described in the text. B. Simulated baroreflex output based on 
input aortic pressure wave obtained from Coleridge et a/. 11 . 

Figure 2B plots model-predicted baroreflex firing rate elicited by 
the normal aortic pressure waveform of Coleridge et al. 12 . The 
model produces the characteristic bursting pattern where peaks 
in baroreceptor afferent firing occur in systole, with the firing rate 
dropping to zero in diastole. 

Responses to ramps of pressure in vivo are compared to data of 
Coleridge et al. in Figure 3. Coleridge et al. adjusted in vivo pressure 
in the aortic sinus by placing hydraulic occlusion at various positions 
along the aorta. Model simulations are based on applying a constant 
increase/decrease to measured P Ao (f), resulting in the pressure time 
series plotted in Figure 3A and 3C. It is assumed that mean pres- 
sure is adjusted at a rate of ±10 mmHg sec 1 , corresponding to the 
experimental measurements. The top panel in the figure corresponds 
to simulations and data associated with the baseline state, where ini- 
tial mean pressure is 100 mmHg. Data on afferent firing rate for this 
experiment are plotted as open circles in Figure 3B. The solid line 
in Figure 3B represent simulated f BR (f) versus P A (t). The red line 
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Table 1 . Model parameters. 


Aorta/large-artery mechanics 




Sensitivity 


V Q = 0.6875 ml 


Aortic mechanics, Equation (1) 


* 


d 0 = 12 mm 


Aortic mechanics, Equation (1) 


* 


C Ao = 0.007 ml mmHg- 1 


Aorta acute capacitance, Equation (2) 


191 


y, 0 = o.4o 


Aorta creep parameter, Equation (3) 


79.6 


T cAo =0 ^ 2SeC 


Aorta creep time constant, Equation (3) 


0.20 


Kinetics of baroreflex afferent firing 


z = 251 .5 sec 

s 


Baroreceptor parameter, Equation (5) 


0.16 


a = 0.0651 sec 1 


Baroreceptor activation rate, Equation (7) 


0.21 


b = 0.2004 sec 1 


Baroreceptor deactivation rate, Equation (7) 


0.35 


S Q = 0.4965 


Baroreceptor saturation constant, Equation (7) 


0.74 


f = 299.8 sec 1 


Baroreceptor gain parameter, Equation (6) 


2.55 


Mechanics of the heart and circulation 


E max = 8 mmHg ml" 1 


Varying elastance heart model, Equation (10) 


* 


E min = 0.25 mmHg ml 1 


Varying elastance heart model, Equation (10) 


* 


T M =0.3 


Varying elastance heart model, Equation (10) 


* 


7^ = 0.15 


Varying elastance heart model, Equation (10) 


* 


H 0 = 75 beat mirr 1 


Varying elastance heart model, Equation (11) 


* 


H, = 100 beat mirr 1 


Varying elastance heart model, Equation (11) 


* 


R out = 1 x10 -4 mmHg min ml -1 


Aortic valve resistance, Equation (12) 


* 


R AQ = 3x10" 4 mmHg min ml" 1 


Aortic resistance, Equation (12) 


* 


R AQ = 0.01160 mmHg min ml" 1 


Large-artery resistance, Equation (12) and Equation (13) 


* 


R v = 3.359x1 0 4 mmHg min ml 1 


Downstream (venous) resistance, Equation (12) 


* 


C A0 = 0.8185 ml mmHg- 1 


Large-artery compliance, Equation (13) 


* 


C w = 329.9 ml mmHg- 1 


Downstream (venous) resistance, Equation (13) 




^ = 625.1 ml 


Unstressed volume of cardiovascular system, Equation (17) 


1.82 


3V=0.40 


Venous creep parameter, Equation (17) 


* 


r . = 120 sec 

cv 


Venous creep time constant, Equation (17) 




a 1 = 0.319 


Arterial and venous compliance parameter, Equation (13) 


0.36 


a 2 = 14.18 


Arterial resistance parameter, Equation (13) 


0.05 


a 3 = 1.036 


Arterial and venous compliance parameter, Equation (13) 


0.40 


a 4 = 2.969 


Arterial resistance parameter, Equation (13) 


0.60 


F 0 = 1125 ml min- 1 


Autoregulation parameter, Equation (14) 


5.96 


F A = 521 .3 ml min- 1 


Autoregulation parameter, Equation (14) 


0.82 


r AR = 6.77 min 


Autoregulation time constant, Equation (14) 


0.20 


r F = 1 5 sec 


Arbitrary time constant, Equation (15) 




Autonomic system 


L =2.76 sec ' 


Baroreflex arc parameter, Equation (18) 


5.5 


Renin-angiotensin system 


t r = 12.61 min 


Time constant for renin production, Equation (19) 


0.30 


r A2 = 1.117 min 


Time constant for angiotensin II production, Equation (20) 


0.30 


t p = 1 5 sec 


Time constant for mean pressure calculation, Equation (19) 




P 1 = 20.12 mmHg 


Steady-state renin-angiotensin system tone, Equation (19) 


0.17 


P 2 = 24.98 mmHg 


Steady-state renin-angiotensin system tone, Equation (19) 


0.37 


g= 245.9 mmHg 


Steady-state renin-angiotensin system tone, Equation (19) 


1.05 


Neurohumoral control of pressure-diuresis/natriuresis 


k A = 0.125 ml sec 1 mmHg -1 


Slope of acute pressure-diuresis relationship, Equation (21) 




P smax = 126.4 mmHg 


Long-term pressure-diuresis relationship, Equation (22) 




P smin = 9.779 mmHg 


Long-term pressure-diuresis relationship, Equation (22) 




0 O = 0.1928 


Long-term pressure-diuresis relationship, Equation (22) 




0 1 = 0.4813 


Long-term pressure-diuresis relationship, Equation (22) 




x k = 10 min 


Time constant for long-term pressure-diuresis, Equation (22) 


0.09 



* parameter not identified based on fitting time-course data; see text for details 
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Figure 3. Baroreflex response to in vivo pressure ramps. A. Applied pressure transients for increasing and decreasing pressure starting 
at a baseline mean pressure of 100 mmHg. Pressure is increased/decreased so that the mean pressure changes at the rate of ±10 mmHg 
sec 1 . B. Simulated baroreceptor firing response to pressure ramps from A is compared to data from Coleridge et alV. C. Applied pressure 
transients for increasing and decreasing pressure starting at a baseline mean pressure of 125 mmHg. D. Simulated baroreceptor firing 
response to pressure ramps from C is compared to data from Coleridge et al?\ See text for details on simulation protocol. In B and 
D, experimental data are plotted as circles; black line represents the simulation prediction; red line represents simulation predictions averaged 
over each heart beat. 

represents f BR (i) and P Ao (t) averaged over each beat. The lower panel 
in Figure 3 corresponds to a hypertensive state where the aortic 
pressure was held at an elevated level with mean of 125 mmHg for 
20 minutes prior to the pressure ramp experiment. Thus the system 
has had 20 minutes to adapt or "reset" to a new mean pressure, 
resulting in a right shift of the baroreflex response/ 5/? versus P Ao . 

The values of the parameters r , 8 Q , f 0 , a, and b were adjusted to 
match the data in Figure 2 A and Figure 3. 

Model Components 3: Mechanics of the heart and 
circulation 

The circulation is modeled as a closed-loop lumped-parameter cir- 
cuit illustrated in Figure 4 A, which ignores the pulmonary circulation 
and treats the heart as time- varying elastance representing the left 
ventricle 15 . The left- ventricular pressure is described by 

P LV (f)=E LV (tyV LV (t), (9) 

where E LV (t) is the left- ventricular elastance and V LV (t) is the volume 
of blood in the ventricle. The elastance is simulated using a smooth 
function that increases in systole and decreases as the heart relaxes: 



E LV (0)= 



((O.75 + 0 w )£ max -£ min ) 



((0.75 + 0 av )£ max -£ min ) 



•(6 



o<e<T M 

+ E imn T M <0<T M + T R , 

T M + T R <e<i 



(10) 



where 6 e (0,1) is the fraction of a heart beat that has elapsed at a 
given time. The parameters E and E . represent minimum effec- 

0 r max mm r 

tive elastance in the ventricle during systole and diastole. The factor 
(0.75 + (j) SN ) multiplying E max accounts for the effect of sympathetic 
tone on contractility, where (j> SN (t) e (0,1) is a model variable (see 
below) representing the sympathetic drive. Under baseline condi- 
tions (j) SN = 0.25 and therefore, under maximal sympathetic stimula- 
tion, cardiac contractility is approximately 175% of baseline. 



The variable 6 is simulated via 



cl6 
dt 



= H, 



(11) 
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Figure 4. Overview of the whole-body model. A. Diagram of cardiovascular circuit model representing the systemic circulation. Flows 
entering and exciting the heart are denoted F jn and F out ; pressure in the aorta and arterial and venous capacitors are denoted P Ao , P A and P v ; 
left-ventricular pressure is denoted P B. Diagram of control systems captured by the model. 



where H = H Q + H x ((j) SN -0.25) is the heart rate, and Q(i) is reset to 
0 each time the value reaches 1. The parameters H Q and H l are set 
to give a baseline heart rate of 75 beats min 1 and a heart rate of 
150 beats min 1 under maximal sympathetic tone. 



where P LV (t) is determined by Equation (9). 

The flows Q and Q . represent the rates of volume uptake/infusion 

*^in ^ urine r r 

and urine production, described below. 



The circuit model is simulated based on equations for the state vari- 
ables 6(t), VJfi, VJtX V A (t), V v (t), VJf), and V sV (t). The governing 
equations for the six volume variables are 



■ max 



= max 



dt 

dV, 
dt 

dV A ^ P Ao ~l 
dt R Ao 

dV v= P A -P, 
dt r AR R A 

dV { 



' P -P ^ 



o, 



p -P 



P ~P 

1 LV 1 Ao 
R out 

p -P 



+ Q r~Q 



dt 

d^y 

dt 



{y s Ao ^sAo)^cA 

={v; v -v sV )/t cV , 



(12) 



The max(0,-) terms in Equation (12) account for the valves, which 
permit flow only in the direction indicated in Figure 4A. 

The neurohumoral control mechanisms captured by the model are 
diagrammed in Figure 4B. The R Ao and R y resistances are set to 
constant values, while other resistances and capacitances vary with 
sympathetic tone and angiotensin II level. Specifically, C A (t), C v (i), 
and R it) are determined by 



c A {t) = 

Cy(t) = 



c A 



(l + a 1 ^ v (O)(l + a 3 0 A2 (O) 
r 



(1 + ^(0X1 + ^2(0) 

^A(0 = ^40 (1 + ®2 W)(l + 0C 4 0 A2 (O), 



(13) 



where C AQ , C yo , and R AQ are constants, and a v a 2 , a 3 , and a 4 are con- 
stants that determine the magnitude of the effects of vasoconstriction 
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via sympathetic tone and angiotensin II. The variable (j) A2 represents 
the plasma angiotensin II activity. The governing equations for (f) A2 
are described below. 

The R A resistance is further governed by a whole-body autoregu- 
lation phenomenon, as incorporated into the Guyton-Coleman 
models 6 ' 8 9 . The function r AR (t) e (0,1) accounts for autoregulatory 
effects on systemic arterial conductivity based on 



Tar dt " M 



1 + tanh 



F - F n 



(14) 



where F(t) is a variable that averages the mean cardiac output with 
a moving average defined by a first-order process: 



dF 

" dt 



r A „R A 



(15) 



The parameters F Q and F l determine the autoregulatory response of 
resistance to changes in flow. The pressures are computed from the 
relationships 



p =(v - V )/r 

1 Ao \ y Ao y sAo)' ^Ao 

p a = v a /c a 

P =(V -V - V )/c 



(16) 



in which venous compliance is simulated using a linear formulation of 
stress relaxation similar to that used for the aorta. Specifically, venous 
stress relaxation kinetics are governed by t cV (V s /dt) = V~ v 



V .„ where 

sV 



v: 



C 



yv ^voi) JVvv ^voi)' 



(17) 



A 120 



105 



90 



75 



60 





f \ Pao I \ 


\.l 


ft [\ 



0.5 



1 

t (sec) 



1.5 



60 



45; 



30 3 

| 

15 




Figure 5. Baseline model operation. A. Model-predicted aortic 
pressure and baroreflex firing rate, obtained with Q jn = 0.5835 ml 
min 1 . This simulation represents a period steady-state of the model, 
in which Q jn = Q out and average pressure is 100 mmHg. B. Model- 
predicted aortic and left-ventricular pressures are plotted for the 
baseline period steady state. 



The constant V yQl in Equation (17) represents an unstressed volume 
for the overall cardiovascular system. 

There are a total of 23 parameters associated with this component 
of the model. Assignment of the values listed in Table 1 was guided 
by a variety of data sets and previous computational models. The 
values of E and E were obtained by scaling the model of 16 to 

max max J ° 

provide reasonable pressure under baseline conditions at ventricular 
volumes appropriate for dog. The cardiac cycle timing parameters 
T M and T R were set to the values used by Beard 16 . Heart rate param- 
eters H Q and H l were determined as described above. The baseline 
resistance and compliance parameters R , R Ao , R AQ , R y , C AQ , C VQ , and 
V VQl were chosen so that under baseline conditions (Figure 5), the 
mean pressure is 100 mmHg, the diastolic and systolic pressures are 
85 and 115 mmHg, and the ejection fraction is 0.58. (The values for 
these seven parameters do not represent a unique set that gives these 
outputs under baseline conditions). The parameters y v and t cV were 
set to match measurement of stress relaxation in the canine jugular 17 . 



The remaining seven parameters in the cardiac and circulation 
model component (a v a v a 3 , a 4 , F Q , F v and t ar ) were identified 
by simulating the responses of the system to volume infusion and 
hemorrhage, as detailed below in Results. 

Model Components 4: Autonomic system 

The whole-body sympathetic tone is represented in the model by 
the variable <j> SN (t) e (0,1) and is determined by the baroreflex arc: 



■*YSN 

dt 



-■fsN^-^-fBR^Sl 



(18) 



Thus, in the absence of baroreflex firing, the sympathetic tone will 
approach the maximum value of 1 . The constant parameter f SN is 
set so that under baseline conditions (j> SN (t) = 0.25. Thus, following 
a rapid severe drop in pressure (j) SN (t) will approach a value that is 
four times the baseline value. 
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Model Components 5: Renin-angiotensin system 

The state of the renin- angiotensin system is captured by two vari- 
ables: <j> R (t) e (0,1) and <j> A2 (f) e (0,1), which represents plasma 
renin and angiotensin II activity governed by the combined action 
of sympathetic tone and pressure on renin release. 

The plasma renin variable (j> R (t) is governed by 



The constant parameters P , P , <b n . and d>. determine how the 
acute pressure-diuresis relationship shifts in response to changes 
in the tone variables (j) SN and (j) A2 . Thus it is assumed that cj> SN and 
(j) A2 have equal and additive effects on renal function. Equation (22) 
assumes that changes in P governed by sympathetic tone and the 
renin-angiotensin system occur with a time constant z k . Further- 
more, the maximal rate of urine production is set to 10 ml miir 1 . 



P dt 



1 - tanh 



PsN F*1 



(19) 



where (j) R decreases with increasing time-averaged arterial pressure 
P. The parameters P l and P 2 determine the shape of the steady state 
relationship between renin production and pressure. This formu- 
lation assumes that the relationship between steady- state (j> R and 
pressure P is shifted by the sympathetic tone. Based on observa- 
tions showing that plasma renin activity and angiotensin II levels 
are nearly perfectly linearly related in vivo 18 , we assumed that <$> A2 {t) 
follows <j> R (i) according to 



(20) 



These expressions simplify the known (and unknown) mechanisms 
governing renin production into a simple phenomenological rela- 
tionship between pressure, sympathetic tone, and the kinetics of the 
renin- angiotensin system. Here the kinetics of the renin- angiotensin 
system are treated with much less detail than other models, which 
account for whole-body pharmacokinetics 19 and saturable enzyme 
kinetics 20 . The five parameters invoked in this model component 
(t r , t A2 , g, P v and P 2 ) are identified based on comparing simulations 
to measurements of pressure, heart rate, and plasma renin activity in 
rabbit during graded hemorrhage, as detailed below in Results. (The 
simulation parameter z p was arbitrarily set at 15 seconds). 

Model Components 6: Neurohumoral control of pressure- 
diuresis/natriuresis 

Regulation of body-fluid volume is assumed governed by a linear 
pressure-diuresis relationship: 



Qur, ne ^ 



min [k^iP-PX lOmlmirT 1 , P>P s 



(21) 



where Q is rate of volume output via the kidneys and k, is the 

urine x J 1 

slope of the relationship between Q . and pressure. The model var- 

x 1 ^ urine L 

iable P s (t) is the variable offset of the pressure-diuresis relationship, 
which is controlled by sympathetic tone and angiotensin II level: 



P -P 

poo p s, max s,min 



1 + tanh 



^SN ' TA2 rp 



dP 



(22) 



The parameters k v P s max , P smin , 0 O> and <j> v were set by matching 
model predictions to data on the pressure diuresis relationship 
under physiological conditions, with angiotensin II infusion, and with 
administration of an angiotensin converting enzyme (ACE) inhibi- 
tor. To compare predictions of Equation (21) and Equation (22) to 
renal output measurements under ACE inhibition, we assumed that 
for the chronic measurements, sympathetic tone was maintained 
at its baseline value and (j) A2 = 0. Fitting data from Hall et al. 21 , it 
is estimated that P~ (<j) SN = 0.25, (j) A2 = 0) = 75 mmHg. Similarly, 
with angiotensin II infusion, we assumed that sympathetic tone was 
maintained at its baseline value and (j) A2 = 1. The data of Hall et al. 
yield an estimate of P~ ((j) SN = 0.25, (j) A2 = 1.0) = 125 mmHg for these 
conditions. Finally, under baseline conditions, with P = 100 mmHg, 
urine output was estimated to be Q = 0.5835 ml min -1 , based 

x urine 

on volume infusion experiments described in the Results. With 
k l = 0.125 ml sec 1 mmHg 1 (see Results), and using the baseline (j) A2 
value of 0.1864 (which is determined from blood withdrawal 
experiments; see Results), it follows that P s °° (<^= 0.25, (j) A2 = 0.1864) 
= 95.3 mmHg. These estimates of P™ at three different values of 
(j) SN + (j) A2 were used to estimate the values of P 



7 P . , <b n , and <f>, . 



The time constant t, was determined from data on the rate of urine 

k 

production following infusion of blood leading to an acute increase 
in pressure and drop in (j> SN and (j) A2 . Analysis of data from these 
experiments is detailed below in the Results together with a direct 
comparison between the data of Hall et al. 21 and model simulations. 

Results 

Results 1 : Response to volume infusion 

Guyton and colleagues conducted experiments in which a large 
amount of blood was infused into an anesthetized dog, resulting 
in a rapid increase in total blood volume of 45% compared to ini- 
tial baseline value. The model of Equation (l)-Equation (22) was 
simulated and compared to data from this experiment to identify 
adjustable model parameters and to probe the predicted response of 
unmeasured model variables to this protocol. 

Figure 6 shows the predicted effects of infusing 45% of initial base- 
line blood volume in a normal animal on arterial pressure, cardiac 
output, rate of urine formation, blood volume, sympathetic tone, and 
renin-angiotensin level together with experimental data obtained 
from Dobbs et al. 22 , Guyton et al. 1 , and Prather et al. 23 on the first 
four variables. (For these experiments Guyton et al. reported a base- 
line mean arterial pressure of approximately 115 mmHg. Since all 
other data sets analyzed here have a baseline pressure of 100 mmHg, 
15 mmHg was subtracted from the data from these experiments to 
achieve consistency). The initial steady state of the model was 
obtained based on a constant infusion of Q in = 0.5835 ml min 1 , which 
matches the mean reported rate of volume output for the baseline data. 
For times 0 < t < 5 min, Q was increased to 0.5835 + 126.82 ml min 1 , 
to result in a total excess volume of 634 ml. After the five-minute 
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infusion, Q was returned to the baseline value of 0.5835 ml miir 1 . 
The model is able to effectively reproduce the trends in the experi- 
mentally measured variables, and predicts that <j> SN and <j> R drop to 
less than one half of their resting values following the volume infu- 
sion in response to the transient increase in pressure. These reduc- 
tions in (j) SN and <j> R confer increases in vascular compliances helping 
to accommodate the substantial volume increase. As volume is 
removed from the system (j) SN and <j> R return toward the baseline val- 
ues, with (j) SN responding substantially faster than <j> R . 



The data plotted in Figure 6, along with data observations on graded 
hemorrhage, were used in identifying several adjustable parameters 
in the model. (See following section for details). 

Results 2: Response to hemorrhage 

Data from Quail et al. 24 on heart rate, pressure, and plasma renin 
activity following graded blood withdrawal in rabbits were used to 
identify parameters associated with the cardiovascular mechanics 
and the renin-angiotensin system. Figure 7 shows data measured in 
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Figure 6. Simulation of arterial pressure, cardiac output, rate of urine formation, blood volume, </> s , v , and ^following infusion of 45% 
of initial baseline blood volume in a normal animal. Data on arterial pressure, cardiac output, rate of urine formation, and blood volume 
are obtained from Dobbs ef a/. 18 , Guyton et al.\ and Prather et al? 9 . The initial steady state of the model was obtained based on a constant 
infusion of Q jn = 0.5835 ml min 1 . For times 0 < t < 5 min, Q jn was set to 0.5835 + 126.82 ml min 1 , to result in a total excess volume of 634 ml. 
The initial condition for the simulation was obtained by setting Q = 0.5835 ml min 1 and running the model to obtain the steady state. 
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Figure 7. Simulation of system response to hemorrhage. The experimental and simulation protocol is to withdraw blood at a rate of 
2% of initial total volume per minute, starting at time 0. The end of the withdrawal period, 17.5 minutes, is indicated in by dashed line in all 
plots. Model simulations are compared to data for mean pressure (A), heart rate (B), and plasma renin activity (D). Panel C plots the model- 
predicted sympathetic tone during the protocol. The initial condition is the same baseline condition used for the simulations of Figure 6. Data 
from heart rate, pressure, and plasma renin activity following graded blood withdrawal are from Quail et a/. 18 with pressure and heart rate 
scaled as described in the text. 



normal rabbits, where the heart rate measurements have been scaled 
to a baseline value of 75 beats miir 1 , renin activity is scaled to 
maximum value of 1, and mean pressure is scaled to an initial 
baseline value of 100 mmHg. In the experiments, blood was with- 
drawn from the animals at a rate of 2% blood volume per minute. 
In model simulations, blood is withdrawn at the experimental rate 
starting at time 0 and with withdrawal stopped at time 17.5 minutes. 
For this experiment, the initial condition was identical to that used 
for the volume infusion experiments of Figure 6, and volume infu- 
sion and urine output were set to zero for the time course. 

Model predictions of H, P, and (j> R compare favorably to experi- 
mental data, with mean arterial pressure dropping from the initial 
value of 100 mmHg to approximately 93 mmHg at the end of the 
blood withdrawal. The observed increase in heart rate is associated 
with an increase in (j) SN from baseline level of 0.25 to 0.81, and an 
increase of (j> R from baseline level of 0.186 to 0.668, at the termi- 
nation of the withdrawal period. During this period, arterial blood 
pressure is protected from a more severe reduction in large part by 
an increase in vascular tone mediated by increases in sympathetic 
tone and angiotensin II level. Simulated cardiac output (not shown) 
also drops during blood withdrawal, reaching a minimum value of 
690 ml mm 1 (approximately 46% of baseline flow) at 17.5 minutes. 

The model is not able to capture the sharp, almost step-like, 
observed increase in renin that occurs following the 10-minute time 



point (corresponding to 20% of blood withdrawn). Specifically, the 
model predicts a more graded response, without the delay observed 
in the data. A better match to the data might be achieved by 
incorporating a delay into the governing equations, or by simulating 
the renin-angiotensin system using a higher-order system of differ- 
ential equations than Equation (19) and Equation (20). The model 
predicts that (j) SN and H peak at the end of the withdrawal period 
and partially recover toward baseline levels within a few minutes 
after the end of blood withdrawal. The renin-angiotensin tone is 
predicted to remain at an elevated level after the withdrawal. The 
partial drops in <j> SN and H are associated with a graded decrease 
in pressure continuing over the entire simulated time course, with 
pressure reaching approximately 85 mmHg at 30 minutes. Pres- 
sure continues to drop after the end of the withdrawal period due 
in part to the action of the whole-body autoregulatory mechanism 
governed by Equation (16) and Equation (17). Via this mechanism, 
the drop in flow associated with hemorrhage causes a reduction in 
arterial resistance, resulting in a recovery in flow and a continuing 
drop in pressure. 

The model predictions in Figure 6 and Figure 7 are most sensitive 
to the values of 13 parameters from the model components of cir- 
culatory mechanics, the renin-angiotensin system, and the kinetics 
of pressure-diuresis: a v a v a y a 4 , F Q , F v t ar , t r , t A2 , g, P v P 2 , and 
t r The values of these parameters were adjusted to match the data 
plotted in Figure 6 and Figure 7. 
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Results 3: Chronic renal function curves 

The steady-state relationship between arterial pressure and rate of 
urine output is obtained from model simulations by varying Q. n , the 
rate of volume infusion, and attaining steady-state model predic- 
tions of steady state pressure where Q =Q . Figure 8 A plots the 
predicted relationship between rate of urine output and mean arte- 
rial pressure under three different conditions: (1) the normal physi- 
ological state; (2) (j> R clamped at 0, representing complete block of 
angiotensin converting enzyme (ACE); and (3) (j> R clamped at 1, 
representing infusion of saturating levels of angiotensin II. Note 
that the pressure-diuresis relationship is plotted in the traditional 
manner with arterial pressure on the abscissa and rate of volume 
excretion on the ordinate, even though in these simulations Q. n is 
varied and steady- state P is computed as a function of Q jn . 

Model predictions are compared to data from Hall et al. 21 on 
steady-state pressure at different levels of salt/volume loading 
in normal dogs, dogs infused with ACE inhibitor, and dogs in- 
fused with angiotensin II. Model simulations effectively match 
experimental data, validating the assignment of parameter val- 
ues described above under "Neurohumoral Control of Pressure- 
Diuresis/Natriuresis". 

Figure 8B illustrates model-predicted steady- state renin and angio- 
tensin II activities (cj> R and (j) A2 ) and sympathetic tone (<j> SN ) as func- 
tions of Q. n for the normal case. (Recall that from Equation (20) 
angiotensin II activity (j) A2 is equal to (j> R in the steady state). As Q in 
is increased (j> R and (j) A2 decrease to maintain pressure at nearly a 
constant level of the simulated range of volume/salt loading. The 
sympathetic tone on the other hand remains nearly constant around 
the baseline level of 0.25. Sympathetic tone remains constant 
because in the model the only determinant of (j) SN is the baroreceptor 
afferent firing rate, which effectively adapts to the small changes in 
pressure that are associated with the simulated range of Q. n . (Thus, 
the model does not capture suppression of sympathetic tone typi- 
cally observed with chronic salt/volume loading). 

Results 4: Response to chronic baroreflex stimulation 

Numerous studies have demonstrated that chronic stimulation of 
the carotid baroreceptor afferent nerve with implantable devices 
results in a prolonged decrease in mean arterial pressure 25 26 . Figure 9 
illustrates a representative data set from Lohmeier et al 21,28 for which 
baroreflex stimulation was applied in dogs for a seven day period, 
followed by several days of recovery. During the stimulation period, 
mean arterial pressure drops by approximately 20%, renin activity by 
approximately 35%, and sympathetic tone (experimentally assayed 
indirectly by plasma norepinephrine) drops by approximately 55%. 

To simulate this experiment, a constant was added to the model- 
predicted baroreceptor firing rate in the equation for sympathetic 
tone. Thus Equation (18) was replaced by 

= fs N (1 {f B R + f«.)K> (23) 

where f stim is a parameter adjusted to match the experimental data. 
Using f stim = 6.5 sec 1 yields predictions that effectively match the 
observations of Lohmeier et al. 21 ' 2 * (Figure 9). Indeed, given that 
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Figure 8. Renal function curves. A. Model predictions are 
compared to data on the steady-state relationship between mean 
arterial pressure and rate of volume infusion (equal to rate of urine 
output) for normal conditions, for angiotensin converting enzyme 
inhibition (0 fl = 0), and for angiotensin II infusion (0 fl = 1). Data are 
obtained from Hall et al. 17 , in which net salt output is reported under 
these three conditions. Rate of urine volume production is assumed 
proportional to rate of sodium excretion, and normalized to the 
rate of urine production at baseline conditions (P= 100 mmHg) for 
the normal case. B. Model prediction for steady-state renin and 
angiotensin II activities (<j> R and (j) A2 ) and sympathetic tone (0 SW ) are 
plotted as functions of Q jn for the normal case (without (j) R clamped). 

only one parameter was adjusted, the comparison of model predic- 
tions to the five variables plotted in the figure represents a strong 
validation test for the model. In particular, the observed reduction in 
renin activity with baroreflex stimulation represents a phenomenon 
that has not been captured by previous modeling efforts 26 28 . Spe- 
cifically, the model of Illiescu and Lohmeier predicts a substantial 
increase in renin production following initiation of stimulation fol- 
lowed by a return to the initial baseline level. In Illiescu-Lohmeier 
model, the transient increase in renin is a result of the drop in pres- 
sure causing an initial increase in renin production. The predicted 
initial increase is later offset by decreased sympathetic nerve activ- 
ity 24 . In our model the observed phenomenon of reduction in plasma 
renin activity emerges as a property of the integrated computational 
model: stimulation of the baroreflex and associated drop in sympa- 
thetic tone is predicted to result in a drop in renin production even 
with the long-term drop in mean pressure. 
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Figure 9. Response to chronic baroreflex stimulation. Electrical 
stimulation of the carotid baroreflex afferent nerve is simulated 
by modifying the normal model by replacing Equation (18) in the 
normal model with Equation (23) during the baroreflex stimulation 
period (for a 1-week period starting on day 0). Data on mean 
pressure, heart rate, urine output, plasma norepinephrine (plotted 
as (j> SN ), and plasma renin activity (plotted as (j) RA ) are obtained from 
Lohmeier et a/. 2324 . 



Thus our model prediction differs from that of the Illiescu-Lohmeier 
model in two critical ways. First, our model captures the observed 
long-term response of renin while the Illiescu-Lohmeier model 
does not. Second, our model predicts no increase in renin during 
the first day of baroreflex stimulation, while the Illiescu-Lohmeier 
model does. Illiescu and Lohmeier predict that an increase in renin 
"was not measured experimentally because initial blood sampling 
was not made until 24 h after initiating baroreflex activation" 24 . In 
contrast, our model predicts that the increase in renin was not meas- 
ured experimentally because it did not happen. 

Results 5: Exploration of the etiology of primary 
hypertension 

Pettersen et al 29 demonstrate that impairment of the baroreflex 
caused by stiffening of arterial vessels represents a viable hypoth- 
esis for the etiology of primary hypertension. Specifically, using a 
closed-loop cardiovascular mechanics model coupled to a detailed 



model of the baroreflex arc 30 , Pettersen et al. 29 show that when the 
stiffness of the aorta is increased to represent changes in mechani- 
cal properties associated with ageing, the model based on their 
mechanogenic hypothesis predicts a substantial increase in mean 
arterial pressure with age. By assuming a constant blood volume 
for all age groups, the model explicitly did not account for the 
regulation of plasma volume and salt through the kidney and the 
renin-angiotensin system following from the hypothesized shift in 
the renal pressure-diuresis/natriuresis function curve. The rationale 
for not including adaptive mechanisms likely to partially ameliorate 
the effects of arterial stiffening on blood pressure was to test the 
explanatory sufficiency of the mechanogenic hypothesis by show- 
ing that it would predict a stronger relation between blood pressure 
increase and age than empirically observed. As a first step towards 
a complete merge of physiological renal function with the model 
developed by Pettersen et al. to expand the prediction space of the 
mechanogenic hypothesis, we studied how the current model respond- 
ed to changes in the aortic compliance. 

Using the current model, which does account for blood volume regu- 
lation by the kidneys, the hypothesis of Pettersen et al. may be fur- 
ther analyzed by determining how simulations of the current model 
respond to changes in the aortic compliance. Figure 10A (solid line) 
plots predicted steady-state mean arterial pressure as a function of rel- 
ative aortic stiffness, C A °/C A , where C\° is the baseline normal value 

7 Ao Ao 1 Ao 

of aortic compliance, and C Ao is the value used to obtain the pressures 
reported in the Figure. The maximum simulated relative stiffness, 
C A °JC Ao = 4, is approximately the average relative stiffness for the 
75-year-old population simulated in Pettersen et al. 29 . This 4-fold 
increase in aortic stiffness increases the predicted mean pressure 
from 100 to 128 mmHg, with systolic/diastolic ratio increasing from 
115/90 to 139/117. Furthermore, over the simulated range of stiff- 
ness model-predicted (j) SN increases from 0.25 to 0.45 and <j> R increases 
from 0.186 to 0.59 (Figure 10B) with a heart rate increase of approxi- 
mately 25% as mean pressure increases from 100 to 128 mmHg. 

As expected, this predicted pressure increase is less than that pre- 
dicted for a model that does not account for blood volume regulation. 
(The model of Pettersen et al. predicts a mean arterial pressure 
150 mmHg for 75-year-old group). The difference between the pre- 
dictions of the current model and that of Pettersen et al. is explained 
primarily by blood volume regulation: with a 4-fold increase in aortic 
stiffness the current model predicts the active blood volume to 
decrease from 1364 to 1050 ml (Figure 10C). This result supports the 
view that normal functioning kidneys are able to partially compensate 
for age-related increases in arterial stiffness 29 . 

Therefore this model demonstrates that it is possible to capture the 
age-related phenotype of primary hypertension as emerging solely 
from changes to the mechanical properties of the aorta and carotid 
arteries. This may seem surprising since the simulation results plot- 
ted in Figure 10 are obtained with mechanical properties of all com- 
ponents of the vasculature other than large vessels associated with 
the baroreflex remaining normal. Indeed, it is expected that vascular 
stiffening would directly impact other physiological processes. For 
example, Beard and Mescam demonstrated that stiffening of renal 
arteries can explain the shift in the acute pressure-diuresis curve 
associated with changes in mean pressure in Dahl S rats 31 . For the 
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Figure 10. Effects of arterial stiffening on mean pressure. A. 

The model predicted for the steady-state mean arterial pressure is 
plotted as a function of relative aortic stiffness, C A °/C A , where C A ° is 

' Ao Ac i Ao 

the baseline normal value of aortic compliance, and C Ao is the value 
used to obtain the simulated pressure. As stiffness is increased (as 
compliance is decreased), predicted mean pressure increases. 
Calculations assume normal salt/volume loading, Q jn = 0.5835 ml 
min~ 1 , resulting in a mean pressure of 100 mmHg at C A °JC = 1. B. 
Model-predicted steady-state sympathetic tone <j> SN , plasma renin 
activity <j> R , and angiotensin-ll activity (j) A2 are plotted as functions 
of C A °JC Ao . C. Predicted blood volume is plotted as a function of 



model results illustrated in Figure 10 (representing normal kidney 
function), the predicted increase in mean pressure associated with 
increased vascular stiffness is partly due to an increase in the pres- 
sure that is required to elicit a given baroreceptor firing rate. As a 
result, the relationship between baroreceptor firing and mean arterial 
pressure is shifted to higher pressure in hypertension, as has been 
observed clinically and in animal models 7 . This phenomenon has been 
called "adaptation" or "resetting". Yet, even though model predic- 
tions are in agreement with the observed adaptation of the baroreflex, 
model predictions are not in agreement with the interpretation that, 
since baroreceptor s adapt, "they cannot participate in the long-term 
control of arterial pressure" 7 . By contrast, the observed adaptation 
is a crucial component of the mechanogenic mechanism simulated 
here. Because stiffening causes a long-term adaptation of the rela- 
tionship between pressure and baroreceptor firing rate, vascular stiff- 
ening can cause a chronic increase in stable mean pressure. In other 
words, since baroreceptor s adapt, they can contribute to the etiology 
of hypertension. 

These model predictions are in agreement with experimental obser- 
vations of Thrasher, who has shown that mechanically suppressing 
baroreceptor sensitivity leads to sustained elevation in arterial blood 
pressure in a canine model 32 33 . Our model predictions are consistent 
with Thrasher's conclusion that mechanical unloading of barorecep- 
tors leads to sustained increases in arterial pressure 34 . 

Parameter sensitivity 

Sensitivity values of parameter estimates were estimated by com- 
puting the relative change in the mean- square error difference 
between data and model simulation associated with a change in the 
parameter value. For example, as explained above, parameters t , 
S 0 , f 0 , a, and b were estimated by matching model predictions to 
data in Figure 2 A and Figure 3. If ^ is the minimum normalized 
mean- square error difference between model prediction and data, 
then the sensitivity associated with parameter r is estimated as |E(r 
+ O.lr ) - E(r ) | /0. 1. Thus a sensitivity value of 1 for a given param- 
eter would mean that a 10% change in the parameter is associated 
with a 10% in the error. Sensitivities for all parameters identified 
based on the data in Figure 1 -Figure 3, Figure 6-Figure 8 are listed 
in Table 1 . Computed sensitivity values range from 0.05 (for a 2 ) to 
greater than 100 (for C Ao ). 

The low sensitivity of results to a 2 can be understood by examining 
the expression of R A (t) in Equation (13). Since the estimated value 
of a 2 is greater than 10, the expression is effectively approximated 
by R A (t) = R AQ a 2 (j) SN (l + « 4 0 A2 ), where the constants R AQ and a 2 cannot 
be separately identified. Since the baseline resistance and compli- 
ance parameters (including R AQ ) are adjusted to (for given values 
of a v a v a v and a 4 ) are adjusted so that the model generates the 
reported baseline pulse pressure, cardiac output, and ejection frac- 
tion, the value of any of a v oc 2 , oc 3 , and/or a 4 becomes arbitrary for 
values significantly greater than 10. 

Summary and conclusions 

A computational model (Source code is permanently available on: 
http://zenodo.org/record/7126 (10.528 l/zenodo.7 126)) was developed 
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and identified to serve as a phenomenological representation of the 
major physiological processes controlling arterial blood pressure. The 
design of the model was based on balancing the compromise between 
complexity/physiological fidelity and identifiability, yielding a model 
of 16 state variables and 46 parameters. Of the 46 parameters, 23 were 
set to previously established values or arbitrary values as explained 
above, and 23 were adjusted to match experimental data illustrated 
in Figure 1-Figure 3, Figure 6-Figure 8. The model captures physi- 
ological phenomena occurring on times scales ranging from millisec- 
onds (e.g. response of baroreceptor firing to arterial pressure) to days 
(e.g. response of mean pressure to chronic stimulation of baroreflex). 
It also effectively captures the response of the cardiovascular system 
to drastic perturbations, including infusion of 45% of blood volume, 
withdrawal of 35% of blood volume, electrical stimulation of the 
baroreflex, and the partial compensatory response of the kidneys to 
blood pressure increase accompanying arterial stiffening. 

The major conclusions from this work are: 

1. The observed adapting/resetting of the baroreflex to long-term 
changes in pressure does not indicate that changes in baroreflex 
function cannot play a role in long-term changes in arterial pres- 
sure. Model simulation reveals how the baroreflex arc and the renin- 
angiotensin system may both contribute in a coordinated manner to 
affect long terms changes in arterial pressure. 

2. It is demonstrated how renal function is not the central or primary 
determinant of long-term arterial pressure in a computer model that 
captures the observed behavior of the renin- angiotensin system, and 
chronic adaptation of the acute pressure-diuresis/natriuresis phenom- 
ena. This result contradicts the classical 'renocentric' view of the 
long-term control of blood pressure even though the model is built on 
the core assumption of the renocentric view that the pressure-diuresis 
relationship acts as a physiological input-output relationship in which 
pressure determines renal output. As a result, it is shown that hyper- 
tension does not necessarily indicate any degree of renal dysfunction. 

3. The physiological response to chronic stimulation of the barore- 
flex can be explained by a mechanism in which arterial blood pressure 
is regulated primarily through the baroreflex arc and the interopera- 
tion of the baroreflex arc and the renin-angiotensin system. 

4. Model simulations are consistent with the mechanogenic hypoth- 
esis 29 that arterial stiffening represents a contributing factor causing 
changes in pressure, sympathetic tone, and renal function associat- 
ed with ageing. Through pressure-diuresis, normal kidney function 
is predicted to be able to partially ameliorate the effects of arterial 
stiffening on blood pressure. 

Regarding the forth conclusion, model simulations predict that stiff- 
ening of the large arterial vessels that are associated with barorecep- 
tors can contribute substantially to long-term changes in pressure. 
This result contradicts the 'renocentric' theory of blood pressure 
control even though the core assumption of the Guyton 'renocen- 
tric' theory of blood pressure control (that the pressure-diuresis 
relationship acts as a physiological input-output relationship in 
which pressure determines renal output) is a core assumption of 



this model. Although it has been argued that because its effective 
gain resets to a given mean arterial pressure the baroreflex cannot 
play an important role in long-term pressure regulation, this analy- 
sis shows how a dynamic system representing the baroreflex and 
capturing the observed resetting phenomenon can play a key role 
in the long-term control of pressure and the etiology of hyperten- 
sion. Model simulations reveal that vascular stiffening increases 
sympathetic tone and shifts the effective baroreflex response to an 
increased pressure baseline. Because these changes cause an increase 
in angiotensin II activity and associated shift in the acute pressure- 
diuresis relationship, the renocentric 7 and mechanogenic 29 hypoth- 
eses for the etiology of primary hypertension could potentially be 
interpreted as compatible. However, as changes to renal function 
arise as downstream consequences of the mechanical remodeling, 
the mechanogenic explanation subsumes more of the biology 
involved in the etiology of hypertension compared with a renocentric 
explanation. This is a major conceptual advance that also provides 
a new interpretational framework for available experimental and 
clinical data. 

The model developed here represents a simplification and several 
components are represented largely in phenomenological terms. 
A more detailed model may be constructed through integration 
of more biophysically based and detailed component models 303135 
into the framework developed here. Regardless, the simplifications 
invoked in the model do not impact the major conclusions. Perhaps 
as important as the simplified modeling framework is the fact that a 
model identified from data obtained from experiments on dogs and 
rabbits, is used to draw insight into human physiology and patho- 
physiology. The implicit assumption here is that arterial pressure 
regulation processes are at least qualitatively similar in dog, rabbit 
and man. 
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Referee Report: 26 November 2013 
SUMMARY 

This work by Beard etaL, is an interesting, consequential and mind-tickling work, definitely 
worth considering and being published, provided that several issues are addressed. 

Strong points: 

• Sound "horizontal" (multi-organ) physiological modelling & simulation (following a 
systems/engineering analysis and conception). 

• Modular, integrative, multi-dynamic approach. 

• Targeted, specific hypothesis - proposes one possible/sufficient etiogenic mechanism in 
the context of age-related hypertension. 

Weak points: 

• Phenomenological equations. 

• Some conclusions drawn are too assertive (especially with respect to such a complex, 
multi-factorial, and evolving, pathology: hypertension); the likeliness of the conclusion 
appears strongly dependent on the limited physiology (and pathology) that is 
implemented. 



1. GENERAL 

In bio-medical research, it has become understood that mathematical modelling and simulation (MS) can 
provide a complementary approach to experiments and are a necessary response to extreme biological 
complexity. Indeed, biochemical and genetic circuits modelling are becoming mature and productive. On 
the other hand, there remains a very serious "modelling gap" to be filled between the ("low-level") 
biological/physiological processes and ("higher-level") pathophysiology. 

Simply put, the "CV-circulatory-pressive" construct by Guyton and colleagues (started in the late 60's, with 
1972 and 1992 milestones; Guyton etaL, 1972; see Montani & Van Vliet, 2009, Moss etaL, 2012) 
remains a landmark, to the point that most, if not all modelling endeavours refer to it in one way or 
another. This is because it was an exhaustive attempt to describe the CV circulatory system in its 
physiological entirety, including high-level interactions such as an integrated CV and renal system. 

Many versions, declinations (and improvements) do exist, including the late HumMod (Hester etal, 201 1 
should perhaps be cited, instead of Abran etaL, 2007, since it is the current evolution of QCP/QHP). 
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However, as mentioned by the Authors, for various reasons it is difficult or impossible to rationally 
understand and exploit those models (e.g. "too large" in the case of HumMod, or "not identified" as in the 
case of Guyton's models of 1 972 or 1 992). 



Beard etal. propose here a "top-down", modular, "engineering" high-quality MS work, as a rationally 
exploitable (numerically identifiable) alternative to Guyton's model, in order to explore hypertension 
etiogenic hypothesis related to arterial stiffness, baroreflex and the renin-angiotensin system (RAS). The 
main incentive relates to the controversed Guytonian view that long-term regulation of arterial pressure is 
driven by the kidney. In that context the Authors try to demonstrate that arterial stiffening (as in ageing) 
and/or the central nervous system (here simplified as the "baroreflex") could provide a sufficient 
explanation(s) of hypertension, through interactions with the kidney pressure-natriuresis function and the 
RAS. 

The significance and relevance of this work for CV regulations and hypertension is obvious, although 

(i) it is limited by the model itself (i.e. effectively implemented features, simplifications, phenomenological 
approach for equations, etc.) and 

(ii) some of the conclusions/assertions should be tempered (among other reasons: given the 
multifactorial/multi-dynamic nature of HT, that it has been resistant to investigation for more than a century 
now, and given the simplicity of the proposed construct). 

2. MAJOR COMMENTS 

In the targeted domain of the pathophysiology of blood pressure regulation and primary hypertension, this 
work proposes a rigorously constructed and parameterized mathematical model of the "long-term control 
of arterial pressure" (but see below). The authors should emphasize however that the model numerical 
instantiation centers on dog (or rabbit). Also of important note (this should be addressed in the revised 
version), no mention is made of the parameterization procedure and tools or sensitivity analysis. 

Moreover, the parameterization is done under the (quite frequent) implicit hypothesis that arterial pressure 
regulation and other CV-related processes only differ quantitatively (e.g. in parameter values) and are 
otherwise qualitatively similar e.g. in dog, rabbit and man. This very point is rarely, if ever, addressed. We 
do feel that in the present turmoil of "physiological" models targeting *human* pathology, this should be, at 
minimum, mentioned and preferably justified. Given that the Authors target human hypertension, and 
especially given the phenomenological approach chosen, why did they not directly parameterize their 
model as a "human" CVR model? 

As a matter of fact, the anatomy and mechanical/physiological/biochemical properties and functions of the 
human CV system, arteries, heart, kidney and RAS are, to the best of our knowledge, now understood 
and numerically characterized enough to allow for such a human quantitative approach to be adopted 
(especially with this articles chosen phenomenological "input-output" approach). This would have given 
some more resonance and relevance to the work. In particular, even if all CV-pressive, renal and 
hormonal processes proved to be the same between all these species, the possibility would remain that 
their relative contributions (and dynamics) to arterial pressure (dys-) regulations are different. 

As their objectives, the Authors aimed to: 

1 . develop a "computer model capturing chronic adaptation of acute renal function in blood pressure 
control" as a more rigorous alternative to Guyton's circulatory model (Guyton etal., 1972) and to 
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Guyton's-derived approaches, particularly SAPHIR (Thomas etai, 2008, - Hernandez etal., 201 1 , 

could be cited too here, as well as Moss etal., 2012) 
2. thus "capturing the relevant mechanisms". 
These are very strong and ambitious objectives, and possibly ab initio flawed by the phenomenological 
(but otherwise legitimate) approach chosen. 

Further, that this is "a model capturing chronic adaptation of acute renal function in blood pressure" 
cannot be claimed with such a simple kidney function, for example: the pressure-natriuresis ignoring 
actions of the ANS (autonomic nervous system) upon renal hemodynamics and tubular function and the 
angiotensin action on glomerular and tubular function (for instance Glomerular Tubulo Balance and 
Na/water reabsorption). Kidney nerve afferents are also ignored. 

What is really addressed in this work are the interactions (possibly etiogenic) between the physiological 
(input/output) modules that are implemented. 

Finally, after more than a century of physiological and clinical research, the etiology of primary 
hypertension remains unresolved; although it is established that (numerous), genetic, environmental, 
age-related and adaptive factors do intervene. Moreover, when in the realm of age-related hypertension 
the Authors do select putative etiogenic processes (arterial stiffening, reduced BR receptors sensitivity) 
and explore potential interactions with RAS and kidney function (phenomenological pressure-natriuresis) 
in a simplified CV construct. 

Driven by the objective (age related HT/arterial stiffening), the proposed model features six central 
physiological functionalities, as modules (model components) implementing: 

MC1 . "Aorta-large artery mechanics" 
MC2. "Kinetics of baroreflex afferent firing" 
MC3. "Mechanics of the heart and circulation" 
MC4. "Autonomic system" 

MC5. "Neuro humoral control of pressure diuresis/natriuresis" 
MC6. "Renin-angiotensin system" 

For each one of these model components (physiological organ/process), the Authors should give, where 
possible, at least one published "first-principle" (physiologically-based) model as a reference, e.g. 
amongst others (i) Moss and Thomas, 2013 for the kidney/pressure-natriuresis under hormonal influence, 
(ii) Guillaud and Hannaert, 201 0 and/or Claassen et al. , 201 3 for the RAS. 

The model is first validated (short-term) vs haemorrhage and then volume expansion. Of note, no real 
validation criterion is given (except for visual adequation through the figures). 

Secondly, the model is validated (long-term = days to weeks) using published steady-state 
pressure-natriuresis kidney function data. Although the Authors target long-term processes, 2 weeks 
(maximal simulated time here) is far too insufficient to refer to the CV, renal and central physiological 
mechanisms truly involved (months-to-years in dogs, years and decades in man). 

OTHER COMMENTS 

• Is one of the input arrows missing in the model figure (fig 4)? 

• A model diagram (rather than the electrical analog presented), as well as for the main hypothesis 
(vs alternative(s)) would be welcome. 
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• To what extent does this work differ from the previous work by the same group (cited, see 
Pettersen etai, 2013)? 

• The Matlab(C) model code and scripts proposed (http://zenodo.Org/record/7126#.UpOLNOJ21 pg) 
have been tested and verified by Dr F. Guillaud (ModTeam/lnserm U1082). The scripts could be 
run, the model could be simulated and results that were visually identical to those in the paper were 
obtained. 

We have read this submission. We believe that we have an appropriate level of expertise to 
confirm that it is of an acceptable scientific standard, however we have significant reservations, 
as outlined above. 

Competing Interests: No competing interests were disclosed. 



1 Comment 



Author Response 

Daniel Beard, University of Michigan, USA 
Posted: 26 Nov 2013 

The reviewers have provided an extensive and insightful commentary on our paper. We agree with 
most of the points raised. We have extracted specific criticisms and responded to them below. 

Weak points: Phenomenological equations; 

This point is addressed in some detail in the responses to the first two reviewers, as well in 
response to specific points below. 

Many versions, declinations (and improvements) do exist, including the late HumMod ( 
Hester etal, 201 1 should perhaps be cited, instead of Abran etal., 2007, since it is the 
current evolution of QCP/QHP). 

The more recent citation has been inserted into the revised text. 

...some of the conclusions/assertions should be tempered... 

The validity of our conclusions is addressed extensively in our response to the first set of 
reviewers. In brief, we agree that some statements needed to be clarified and rephrased. See 
responses above for details. 

The authors should emphasize however that the model numerical instantiation centers on 
dog (or rabbit); the parameterization is done under the (quite frequent) implicit hypothesis 
that arterial pressure regulation and other CV-related processes only differ quantitatively 
(e.g. in parameter values) and are otherwise qualitatively similar e.g. in dog, rabbit and 
man. 

This is an important point that cuts across all areas of computational physiology. Large-scale 
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models are built based on data from multiple sources, multiple species, different conditions (in 
vitro, in vivo, etc.) In the end, all of our models are chimeras. We have taken the liberty of adapting 
the reviewers' language at the end of the discussion section of the revised paper. (We do point out 
that our implicit assumption is the central assumption of all biomedical research using animal 
models. It might not be a perfect assumption, but it does give us a starting point.) 

Also of important note... no mention is made of the parameterization procedure and tools 
or sensitivity analysis. 

Please see our detailed response to this point in the response to the first review. 

For each one of these model components (physiological organ/process), the Authors 
should give, where possible, at least one published "first-principle" 
(physiologically-based) model as a reference 

This is an excellent idea in principle. However, it is not clear if it is possible. For example, we are 
unclear as to how the Moss and Thomas paper could be characterized as representing the 
hormonal influence on pressure-natriuresis from first principles. Specifically, the Moss and Thomas 
paper represents the action of angiotensin-ll "via the parameter FrANG, which adjusts the 
response of the pressure-natriuresis mechanism." The parameter FrANG appears in the model 
equations as a factor multiplying a phenomenological expression for the fractional reabsorption, 
and the effect of the hormone is simulated by varying FrANG between "a maximal value of 1 
(minimal angiotensin II) and a minimal value of 1/6 (maximal angiotensin II)". Thus, the model is 
every bit as "phenomenological" as the one presented here. (This point is not to be taken as a 
criticism of the excellent study of Moss and Thomas.) 

The papers by Guillaud and Hannaert, and Claassen et al. do account for the RAS system with 
more physiological fidelity/detail than our model. We do disagree with the assertion that any of the 
models discussed here are built from "first" principles. All of these papers invoke the same sort of 
phenomenological expressions for certain components. 

Although the Authors target long-term processes, 2 weeks (maximal simulated time here) 
is far too insufficient to refer to the CV, renal and central physiological mechanisms truly 
involved... 

The reviewers may be correct that our model is far too insufficient to provide insight into a 
pathophysiological mechanism that emerges over decades. Yet, we refer the reviewers (and the 
readers) to the stated conclusions in the paper. 

Is one of the input arrows missing in the model figure (fig 4)? 

The net volume output rate is Q_urine - Q_in. The figure has been revised to clarify this. 

A model diagram (rather than the electrical analog presented), as well as for the main 
hypothesis (vs alternative(s)) would be welcome. 

This has been added. (See Figure 4B.) The main hypothesis here is that pressure-diuresis acts as 
a physiological input-out relationship, and that renal function is controlled through the RAS and the 
autonomic system. Thus, fundamentally, the model is constructed to capture the Guyton concept 
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of the long-term control of arterial pressure. (Therefore, it is not clear what alternatives we might 
indicate in a model diagram.) 

To what extent does this work differ from the previous work by the same group..? 

Please see the Section "Results 5: Exploration of the etiology of primary hypertension" in the 
revised text. 

Competing Interests: No competing interests were disclosed. 




Johnny Ottesen 

Department of Science, Systems, and Models, Roskilde University, Roskilde County, Denmark 



Approved: 14 November 2013 
Referee Report: 14 November 2013 

Complex physiologically systems are often subdivided into parts and studied separately by experimental 
methods. But in experimental physiology, subsystems are not easily put back together and conclusions 
made concerning subsystems have limited validity for the original undivided system as their function and 
dynamics may have changed when isolated from one another. Parts that cannot be isolated 
experimentally may be studied (separately) using modeling. Mathematical modeling is able to unfold the 
influence that each of the processes has on the overall dynamic and thus modeling is highly adept at 
separating systems into components simple enough for their structures and functions to be studied in 
isolation. Mathematical modeling is the only controlled way to put the pieces back together by using 
equations that represent the system's components and processes, as well as its structures and 
interactions. 

In "A computational analysis of the long-term regulation of arterial pressure" the authors follow this 
strategy nicely in a successful attempt to investigate what they denote the c renocentric' theory by Guyton. 
Here the authors combine cardiovascular and baroreflex models with models of the renin-angiotensin 
system whereby some important conclusions are reached. Among these conclusions the authors 
emphasize that the baroreflex may play a key role in the long-term regulation of the arterial pressure and 
especially in hypertension in contrast to the Guytonian view. 

The authors follow the physiologically structure well in their models as they carefully validate the different 
sub-models/elements by comparing them with existing various animal experiments. These sub-models 
are then integrated into a well-behaved and solid overall model. 

However, there are some methodic drawbacks I would like to emphasize: Most model equations are in 
some sense ad-hoc choices (or phenomenological, as the authors mention at the very end). These 
equations are not based on first principles and hence they do not represent the underlying mechanisms. 
Instead the approach is to assume a first order transient toward some static sigmoidal relationship. In 
addition, the authors make a lot of self-references, which shadows the fact that many of the used 
sub-models are adjustment taken from elsewhere and ideas from other sources. Strengthening the paper 
with such information would definitely support the approach taken. 

The paper is very well structured and written but a few misprints have sneaked in: 
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Formula (3): Should the fraction be / C A instead of C A /? 

N 7 Ao Ao 

Formula (4): An index o is missing on V gA 
On page 4, coulomb 2, line 2: should be 

On page 4, coulomb 2, line 5: the population transition from an active state at a rate . . . ' Please state 
reference (or state that this is a hypothetical assumption)! 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 

Competing Interests: No competing interests were disclosed. 



1 Comment 



Author Response 

Daniel Beard, University of Michigan, USA 
Posted: 26 Nov 2013 

Thank you for the attention to our paper and the constructive criticisms. 

Most model equations are in some sense ad-hoc choices (or phenomenological, as the 
authors mention at the very end). These equations are not based on first principles and 
hence they do not represent the underlying mechanisms... In addition, the authors make a 
lot of self-references... 

We agree with the reviewer's assessment. To reiterate a statement from the paper summary, the 
model represents a balance complexity/physiological fidelity and identifiability. The sensitivity 
analysis added to the revised paper indirectly addresses the issue of the phenomenological nature 
of model components. Finally, we have revised the text to clearly state the nature of the model 
earlier than the last paragraph of the paper. See third paragraph on page two, section on "Model 
Components 5", and first paragraph of the summary. 

Also, we point out that while we have cited some previous modelling works, the majority of 
components used here have been constructed de novo, with the overall goal in mind. The varying 
elastance heart model is adapted from the cited reference. And of course there is nothing new 
about lumped-parameter models of the circulation. But most of the other (phenomenological) 
formulae are constructed to match the observed dynamics of responses of these systems to the 
experiments used to identify the model. 

Formula (3): Should the fraction be / C. instead of C. /? 



No, it is correct in the paper. One way to see this is to think about the limit where C Ao goes to 

infinity. In this case the vessel becomes infinitely compliant if we wait long enough and thus does 
not support any finite pressure. 
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Fixed. 

On page 4, column 2, line 2: should be d Q 

Fixed 

On page 4, column 2, line 5: the population transition from an active state at a rate 
Please state reference (or state that this is a hypothetical assumption)! 

The assumption of Equation (7) is stated: "It is assumed that the baroreceptors within the 
population transition from an active to inactive state at a rate proportional to the firing rate". We do 
not know what factors are responsible for the slow "resetting" of the baroreflex response. But, it is 
compelling that this simple equation is able to capture the nonlinear phenomena of Figures 2 and 
3. 

Competing Interests: No competing interests were disclosed. 



John Osborn 1 , Viktoria Averina 2 

1 Department of Integrative Biology and Physiology, University of Minnesota, Minneapolis, MN, USA 

2 University of Minnesota, Minneapolis, MN, USA 

Approved with reservations: 04 November 2013 
Referee Report: 04 November 2013 

This paper by Beard etal makes three major claims: 1) to provide a mathematical description for Guyton's 
existing model; 2) to show how stimulation of the baroreflex arc can explain long-term changes in blood 
pressure via the renal function curve employed in Guyton's model; and 3) to show how arterial stiffening 
can lead to hypertension without renal injury. I commend the authors for making a much needed 
mathematical contribution to the field of computational physiology, where one of the most widely 
recognized and used models of long-term control of arterial pressure has only very rudimentary 
mathematical description and analysis published. Moreover, the authors are bridging the gap between the 
renocentric view of hypertension etiology and the emerging experimental data that points to an 
independent long-term role of neural regulation. For these reasons the paper is a valuable addition to the 
existing body of knowledge. 

There are several issues however, that need to be addressed in a revised version of the paper. The 
overarching recommendation is to correct the paper's focus and claims to make them more consistent 
with the model design, as well as to make the model's description easier to understand. In particular, I 
believe the model should not be presented as a generic model with the "representation of the major 
physiological processes controlling arterial blood pressure" but instead be positioned as a conceptual 
framework which aims to illustrate how the baroreflex arc may be involved in long-term blood pressure 
control. Phenomenological mathematical representation, a small amount of detail, limited analysis, and 
quality of shown simulations are appropriate for a conceptual model but are insufficient for a model from 
which one could be extracting conclusions on the particular physiological mechanisms to be implicated in 
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hypertension development. 
Suggested revisions: 

There are several versions of Guyton's model ranging from very simplistic to very complex 1 . Some 
models are extremely detailed but have a closed architecture (e.g. hummod.org). The model created by 
the authors is not sophisticated enough to correspond to most of these versions. Please state clearly 
which version of Guyton's model is being considered here. In my opinion, the model proposed here 
seems of approximately the same level of complexity to the one presented by Guyton and Coleman 
'Long-term regulation of the circulation: interrelationships with body fluid volumes' in the text Physical 

bases of circulatory transport: Regulation and exchange. (1 967) 2 

1 . Focus the goals according to what was accomplished in the paper: 

• "capturing the relevant mechanisms in an identifiable mathematical model" is too 
over-reaching since the three key equations (18, 19, and 21) do not capture mechanisms 
but are of phenomenological nature. 

• "probing hypotheses regarding the long-term control of arterial pressure" - the model is both 
too limited and phenomenological to allow for any testing of various hypotheses. It is, 
however, fit for providing a possible explanation for two hypotheses: roles of the baroreflex 
involvement and arterial stiffening. Otherwise, for testing, one would require a more 
mechanistic description, inclusion of other key organs and regulatory mechanisms, and a 
more rigorous mathematical analyses. 

2. The authors provided an explanation of how the baroreflex may be involved in long-term control of 
blood pressure. This is a novel contribution to the field and should be the focus of the discussion: 

• Add the discussion and references of experimental knowledge of baroreflex stimulation, 
inhibition, and especially long-term resetting. 

• Describe how Guyton's models incorporate neural control and explain the above finding, 
especially the effects of chronic stimulation and how it is different from the proposed model 

(e.g. 3 ) 

• Discuss existing knowledge or a possible physiological explanation for the finding that renin 
activity will be diminished with baroreflex stimulation. 

• Discuss what other important physiological mechanisms are not reflected in the model and 
how they may play a role. In particular, explain why macula densa sodium delivery is 
thought not to be important enough to include in the model, despite the fact that it is known 
to be one of the key regulators of renin activity. 

• Do not confuse the model's predictions with its assumptions when discussing simulation 
results. The model assumes and is explicitly designed to tie the baroreflex, renin, and SNA 
with long-term pressure control. An explicitly designed assumption is not the same as 
prediction. Thus, statements like Conclusion 1 in "Summary and Conclusions" are 
inappropriate. The fact that the model's assumptions allow for physiologically plausible 
simulations merely supports the notion that baroreflex may be important in the way the 
authors have described it. 

3. Improve the clarity of the mathematical description: 

• Provide motivation for the chosen elements of the model and outline the model in its entirety 
first. For example, a flowchart or a diagram of the model may help readers to grasp the 
overall idea. 

• Add at least a basic analysis of parameter sensitivity, with clear description in the Methods 
section. 
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• In Table 1 please provide physiologically-based description of parameters rather than 
referring to the equations where they are used (e.g. Emax is the maximal left-ventricular 
elastance, rather than "Varying elastance heart model, Equation (10)"). Otherwise, the 
definitions are circular: parameters are used in the equations without an explanation of their 
meaning in either the text or the table. 

References 

1. Guyton AC, Coleman TG, Granger HJ: Circulation: overall regu\a\\on.Annual Review of Physiology. 
1972; 34: 13-44 PubMed Abstract I Publisher Full Text 

2. Guyton AC, Coleman TG, Granger: Lowering of blood pressure during chronic suppression of central 
sympathetic outflow: insight from computer simulations. Physical bases of circulatory transport: Regulation 
and exchange. 1976. 

3. Illiescu R, Lohmeier TE: Lowering of blood pressure during chronic suppression of central sympathetic 
outflow: insight from computer simulations.C//r?/ca/ and Experimental Pharmacology and Physiology. 
2010; 37 (2). PubMed Abstract I Publisher Full Text 

We have read this submission. We believe that we have an appropriate level of expertise to 
confirm that it is of an acceptable scientific standard, however we have significant reservations, 
as outlined above. 

Competing Interests: No competing interests were disclosed. 



1 Comment 



Author Response 

Daniel Beard, University of Michigan, USA 
Posted: 26 Nov 2013 

We are grateful to the reviewers for paying such serious attention to our paper and offering many 
constructive criticisms. While we have revised the paper to address the concerns and clarify the 
presentation, we disagree with the general interpretation of our model as "a conceptual framework 
which aims to illustrate how the baroreflex arc may be involved in long-term blood pressure control. 
" As stated in the summary, the model "was developed and identified to serve as a representation 
of the major physiological processes controlling arterial blood pressure." For example, there is no 
equation or parameter value artificially injected into the model to require pressure-diuresis to not be 
able to return the pressure to a baseline value during baroreflex stimulation. In fact, we could 
parameterize the model so that the chronic pressure-diuresis phenomenon acts as an 
"infinite-gain" control system, always returning pressure to its set point. But in doing so, the model 
would not be able to match the data used to identify it, or to predict the data used to validate it. 

Regarding the level of detail captured by the model, the reviewers are correct that at some level, all 
components of the model are phenomenological in nature. (This is true of all mathematical 
models.) Yet regarding the question around the level of detail and complexity of the model, we 
stand by the assertion that this model represents the most detailed model of long-term control of 
arterial pressure that is identified in the literature. (Here we are invoking the word "identified" to 
indicate that the model equations and associated parameter values are reported along with 
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comparisons of model simulations of data used to identify the parameters.) 



In fact, one of the unstated goals of this paper is to inject some clarity and rigor into this field with 
precise and concrete language. Therefore we do not know how to respond to the statement that 
our model "is not sophisticated enough to correspond to most of the versions of the Guyton model 
because it is not clear, precisely, what is meant by "sophisticated" or "correspond". In discussing 
the Guyton models, we clearly state that: "The 1972 realization of the model invokes approximately 
160 variables, and several hundred adjustable parameters. More recent versions of the model 
involve many thousands of variables and tens of thousands of parameters ." Although it is true that 
our model has less variables and parameters than any version of the Guyton model that we are 
aware of, is the sophistication of a model measured purely by the number of variables or adjustable 
parameters? Since no version of the Guyton model captures pulsatile mechanics, no version of the 
model can capture the baroreceptor resetting phenomena illustrated in Figures 2 and 3 of our 
paper. Nor does the HumMod model capture the system response to baroreflex stimulation 
illustrated in Figure 9. Therefore one could argue that our 16-variable model is more sophisticated 
than the 5,000-variable Guyton model. 

Regarding the specific suggestions: 

"capturing the relevant mechanisms in an identifiable mathematical model" is too 
over-reaching since the three key equations (18, 19, and 21) do not capture mechanisms 
but are of phenomenological nature. 

More specifically, we stated that one of the "goal" of the study is "capturing the relevant 
mechanisms in an identifiable mathematical model ." We believe we have done that. However, in 
the revised paper we have changed the goal to "representing the relevant mechanisms in an 
identifiable mathematical model" 

"probing hypotheses regarding the long-term control of arterial pressure" - the model is 
both too limited and phenomenological to allow for any testing of various hypotheses. 

Again, we counter that we have met the stated goal. We are not claiming to have proved any 
hypothesis. Yet we have introduced some important disproofs, as outlined in some of our major 
conclusions. Specifically, it is stated in the literature as 'scientific facts' that; the renal function 
curve represents the dominant long-term controller of blood pressure and, that since the baroreflex 
"resets" it cannot play a role in chronic changes in mean pressure. Our model analysis has cast 
significant doubt upon the validity of both of these hypotheses. Admittedly, these findings 
represent somewhat of a double negative i.e., a model capturing baroreflex resetting is used to 
demonstrate how a baroreflex that resets could play a role in long-term changes in pressure. 
Therefore, we have shown that the fact that the baroflex resets, does not mean that the baroreflex 
cannot play a role in long-term pressure control. Therefore we are confident that it is appropriate to 
state that we have constructed a model with the goal of "probing hypotheses regarding the 
long-term control of arterial pressure." 

Add the discussion and references of experimental knowledge of baroreflex stimulation, 
inhibition, and especially long-term resetting. 

Data in Figures 2 and 3 represent direct experimental interrogation of baroreflex stimulation, 
inhibition, and long-term resetting. The resetting phenomenon is also indirectly probed in Figure 6 
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where, over the time course of the experiment, the baroreflex response to the acute pressure 
increase gradually diminishes. The phenomenon and experimental observations are discussed in 
association with our simulations of these experiments. These are also discussed in the section on 
arterial stiffening and the etiology of hypertension. Finally, we have added a couple of sentences 
(end of Results Section 5) comparing our model predictions with the long-term observations of 
Thrasher on baroreceptor unloading. 

Describe how Guyton's models incorporate neural control and explain the above finding, 
especially the effects of chronic stimulation and how it is different from the proposed 
model... 

A discussion on this point has been added to the section "Results 4". 

Discuss existing knowledge or a possible physiological explanation for the finding that 
renin activity will be diminished with baroreflex stimulation. 

It is not our finding that renin activity will diminish with baroreflex stimulation. That is the prior 
experimental observation. The mechanism in the model is through the observed dependence of 
renin release on sympathetic tone. At a deeper level, the mechanism may simply be explained as 
follows: decreased sympathetic tone leads to decreased renal arterial tone, leading to increased 
glomerular filtration and/or medullary blood flow, leading to increased flux of sodium to the macula 
densa, leading to decreased renin production. While we would rather leave such speculation out of 
the manuscript, it is published here as part of the open review process. 

Discuss what other important physiological mechanisms are not reflected in the model 
and how they may play a role... explain why macula densa sodium delivery is thought not 
to be important enough to include in the model, despite the fact that it is known to be one 
of the key regulators of renin activity. 

It is not that we believe that "macula densa sodium delivery is. ..not. .important enough to include in 
the model." Rather the phenomenon is captured in a simplified way. To clarify this point, we have 
added the following text to the section on "Neurohumoral control of pressure-diuresis/natriuresis": 

"These expressions simplify the known (and unknown) mechanisms governing renin production 
into a simple phenomenological relationship between pressure, sympathetic tone, and the kinetics 
of the renin-angiotensin system." 

A more biophysically detailed model would be required to represent the direct link between sodium 
transport and renin production. At the bare minimum, such a model would represent the 
biophysical processes linking pressure and solute transport in the kidney, including effects of 
sympathetic innervation on arterial smooth muscle stress development, mechanics of blood and 
filtrate flow, solute advection, diffusion, permeation, exchange, and pumping. 

Do not confuse the model's predictions with its assumptions when discussing simulation 
results... Conclusion 1 in "Summary and Conclusions" is inappropriate. The fact that the 
model's assumptions allow for physiologically plausible simulations merely supports the 
notion that baroreflex may be important in the way the authors have described it. 

We agree that this first stated conclusion does not represent a conclusion and what we are trying 
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to say is now better captured in the revised text for Conclusion 1 . 



This same sentiment applies to the other conclusions, which are theoretical. For example, we 
conclude that we have demonstrated that "that hypertension does not necessarily indicate any 
degree of renal dysfunction." It may seem like an equivocation to insert the word "necessarily". 
However, the literature is full of claims that due to the phenomenon of pressure-diuresis, 
hypertension necessarily does indicate some degree of renal dysfunction. Our theoretical analysis 
however would disprove that widespread claim. 

Provide motivation for the chosen elements of the model and outline the model in its 
entirety first. For example, a flowchart or a diagram of the model may help readers to 
grasp the overall idea. 

We have added a diagram illustrating the control processes captured by the model. (Figure 4B.) 

Add at least a basic analysis of parameter sensitivity, with clear description in the 
Methods section. 

This is an excellent idea. A detailed sensitivity analysis is not feasible because the model is 
expensive to simulate (particularly for the long time frames of Figure 9, for example), the number of 
parameters is large and would require calculation of 46x46 Fischer Information matrix (or 
cross-correlation matrix), and we do not have standard error estimates for all of the data used to 
identify the model. However, as suggested we have done a basic analysis, and reported a 
sensitivity index in Table 1 . Results are discussed in a new section on parameter sensitivity. 

(Also, the sensitivity analysis facilitated as a slightly better parameterization of the model. The 
results have change imperceptibly.) 

In Table 1 please provide physiologically-based description of parameters rather than 
referring to the equations where they are used... Otherwise, the definitions are circular: 
parameters are used in the equations without an explanation of their meaning in either the 
text or the table. 

We believe that it useful to the reader to indicate in the table where the parameter is used in the 
model. Therefore, the indication of which equations given parameters appear in has been retained 
in the table. In addition, descriptions of the meaning of all parameters have been added wherever 
they were missing in the text. 

Competing Interests: No competing interests were disclosed. 
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